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1 Introduction 



Recent studies of type la Supernovae P have demonstrated that the expansion of the 
universe is accelerated. In the context of Einstein gravity this can only be due to an 
exotic form of energy called dark energy or quintessence with negative pressure, which 
dominates the cosmic dynamics The best fit of the observational data from cosmic 
microwave background experiments [3^ fixes the energy density of dark energy to be 
Qde = 0.7 in units of the critical density, with matter accounting for Qm — 0.3, mostly 
composed of dark matter. The best fit of the observational data also favours a time- 
varying form of dark energy |3] and there is marginal evidence for a negative equation 
of state parameter w = P/ p < —1, where P and p are the energy density and pressure 
of the dark energy, respectively |4j. 

Most models of dark energy are based on scalar fields, which are natural candidates 
from the theoretical point of view because it is well known from the inflationary theory 
of the primordial universe that a scalar field slowly rolling in a nearly fiat section 
of its self-interaction potential V (0) can fuel a period of accelerated expansion of the 
universe. However, a single canonical scalar field in Einstein gravity cannot explain the 
range of values w < —1 of the equation of state parameter preferred by the observations 
EI- In fact, a canonical scalar field in general relativity has energy density and 
pressure 

P=y + m, P=^-n0), (1-1) 

where an overdot denotes differentiation with respect to the comoving time t. In the 
spatially fiat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

ds"^ = -dt^ + a^{t) {dx^ + dy"^ + dz^) (1.2) 

in comoving coordinates {t,x,y,z), the Einstein- Friedmann equation 

H = -^{p + P) (1.3) 

(where k, = 8tt G and H = a/ a is the Hubble parameter) implies that w < —1 is 
equivalent to if > 0, a regime called superacceleration as opposed to acceleration d/a = 
H + H"^ > Then, eq. JOl) yields H = -k4>'^/2 < 0, with the upper bound of the 
inequality being attained by de Sitter spaces with H =const. A regime with H > 
(or w < —1) requires an extension of the theory [H]. One possibility that has received 
much attention in the literature jHl El ^| is that of phantom fields, i.e., the scalar 
is formally allowed to have a negative kinetic energy and 0^/2 is replaced by —0^/2 in 
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eq. (|1.3p . This modification allows for dynamics capable of breaking the w < —1 barrier. 
Motivations for phantom fields arise in string/M theory and in supergravity [llj. 

Another possibility |T2] is to allow the scalar (p to couple nonminimally to the Ricci 
curvature through the action 



S 



(fx 







;i-4) 



Both phantom and non-minimally coupled scalar field models can be seen as scalar- 
tensor theories of gravity. Scalar-tensor models of dark energy (called extended quin- 
tessence) were considered even before the observational evidence for superacceleration 
was reported jT2] (see Refs. [ini El for reviews of cosmology in scalar-tensor gravity). 
Perhaps the most striking consequence of superacceleration is that the universe may end 
its existence in a Big Rip singularity at a finite time in the future p. Whether a Big 
Rip is unavoidable or another fate awaits the late universe is determined by the presence 
and size of the attraction basins of Big Rip solutions and/or other late time attractors 
in the phase space of the dark energy model. In certain models the universe eventually 
exits the superacceleration regime and converges to a de Sitter attractor, expanding for 
an infinite time [12]. In other models the Big Rip is the only late time attractor (e.g., 

da nil). 

In this paper we study the stability of de Sitter solutions (when they exist) in ex- 
tended quintessence. The generalized gravity action 



S 



f{^,R) ^(0) 
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contains scalar-tensor gravity as the special case in which the function / (0, R) is linear 
in R. The field equations in the FLRW metric p.2|) are 
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H = -—{uj<P' + F-HF 



;i.6) 
;i.7) 
;i.8) 



where F = df/dR. The Hubble parameter H is a cosmological observable and the 
scalar is a natural dynamical variable to consider, hence we take H and as physical 
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variables. Then the only equilibrium points of the dynamical system (jl.6|) - (|1.8|) are 
de Sitter spaces with constant scalar field {Hq, 0o) and they exist subject to the two 
conditions 

6H^Fo-fo + 2Vo = 0, (1.9) 



fo' - 2K)' = , (1.10) 

dV 



and 

((/-o.^Jo)' 



where F^^F (Ho, 0o), fo^f (0o, Ro), Vo ^ V (0o), = f , fo' ^ % 

4)0 

/ i = 12Hq. Equations ()1.9p and (jl.lOj) are straightforward to derive from eqs. p.6p 
(II. 8p and appeared in the literature before for various special choices of the function 
/ (0) R) [IH|-|ni]- There are only two conditions because only two equations in this set 
are independent. The stability of these de Sitter spaces is usually assessed with respect to 
homogeneous (time-dependent only) perturbations; in Ref. [221 we studied stability with 
respect to inhomogeneous (space- and time- dependent) perturbations, to linear order. 
The difficulty with inhomogeneous perturbations is related to their gauge- dependence 
and to deal with this problem we used the covariant and gauge-invariant formalism of 
Bardeen, Ellis, Bruni, Hwang and Vishniac in the formulation given by Hwang and 
Hwang and Noh for generalized gravity We found that, when 0o 7^ 0, the de Sitter 



fixed point {Hq, 0o) = ^"qf^^ 'i 0oj is stable if and only if 



— ^0 ^ F\, 



-0 (l + St) 



<0 (1.11) 



(see Refs. jl8 | IT? H 1 ^ 1^ for similar conditions in generalized gravity). A different for- 
mula applies if 00 = but we will not be concerned with this case here, while contracting 

de Sitter fixed points with Hq = ~\J ^^^p^ are always unstable [22] • We refer the reader 
to the Appendix for a sketch of the derivation of the stability criterion (|1.11|) . and to 
Ref. for a detailed discussion. 

Our goal is to apply the stability condition p.ll|) to extended quintessence scenarios 
that appeared recently in the literature and to scalar-tensor infiationary scenarios of the 
early universe j2l 
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2 Non- minimally coupled quintessence 



We begin with the extended quintessence model of Carvalho & Saa based on a 
non-minimally coupled scalar field and described by 



/(0,/?) = (i-e0')i? 



UJ 



V{(f)) =Ae- 



(2.1) 



where ^ < is a dimensionless coupling constant, A and a are positive constants, and 
units are used in which k = Stc G = 1. Equations ()1.9p and p.lOj) for the existence of 
de Sitter fixed points {Hq, (po) yield 



aAe-"'' 



3(1 12^00 




(2.2) 



which correspond to eqs. (14) and (13) of Ref. j27|, respectively. These authors consider 
the range of values of the coupling constant ^ < — (T^/4, and 0^ < 0. The criterion for 
stability ()1.11|) in the theory ()2.1|) reduces to the one previously found in Ref. j2Hl and 
can be written as 



00 



X = ^00 



l-3x 
1 — a; 



(2.3) 



if 00 7^ 0. It is convenient to use the variable y = ^ with < y < 1. In the case of the 
fixed point 



(-))2 
) 



3 [1 - e(0 

the stability condition ()2.3|) simply reduces to 



a 




(2.4) 



y + 



y>o 



(2.5) 



which is always satisfied, hence this equilibrium point is stable against linear inhomo- 
geneous perturbations. This conclusion complements the work of Carvalho and Saa 
fI7\ who consider only homogeneous perturbations but to arbitrary order by exhibiting 
a suitable Ljapunov function. On the one hand their approach provides information 
on the attraction basin; on the other hand, our approach takes into account the more 
general inhomogeneous perturbations by tackling the gauge-dependence problems. Re- 
garding the parameter space, Carvalho and Saa only prove stability analytically for 
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—5/6 < ^ < — cr^/4, although they conjecture stabihty for a much larger range of 
values of this parameter on the basis of numerical integrations probing down to val- 
ues ^ < —100. Our approach instead establishes linear stability analytically for any 
^ < — cT^/4. It is to be noted however that the values of a and ^ are not arbitrary: 
they are subject to bounds placed by observational constraints on the variability of the 
effective gravitational coupling Gejj(0) = G (1 — Svr 0^) ^ [2^]. 

The second fixed point ^-f^o? ^o^'' j is found to be unstable against homogeneous 
perturbations in Ref. [27]; this is of course true also for inhomogeneous perturbations. 
In fact the stability condition ()2.3p reduces to y{l — y) < 0, which is never satisfied 
because y only assumes values in the range < y < 1. 

The constraints on the variation of Geff mentioned above [29j do not apply to in- 
flation but they are relevant to more recent epochs of the cosmic history (e.g., during 
galaxy formation) and to scenarios of dark energy in the present universe. Although 
most constraints depend on specific quintessence models or parametrizations of the vari- 
ation of Geffi the general feeling is that it is rather difficult to reconcile the available 
constraints on the variation of Geff with a successful extended quintessence scenario. 
While the latter requires low energy scales, high energy scales are generally required to 
implement the modest variations of G^ff allowed. The realization of a successful ex- 
tended quintessence scenario compatible with all the observational constraints and free 
of fine-tuning problems may in the end reveal to be very problematic [22] and this issue 
deserves further investigation. 



3 Inflation with non-minimal coupling 

The approach presented here is suitable to study also the stability of de Sitter-like 
infiation in scalar-tensor gravity; in this context the presence of a de Sitter attractor 
justifies the use of the slow-roll approximation to infiation in the early universe We 
turn our attention to the recent non-minimally coupled scalar field infiationary models 
of Refs. (see also 1331) which employ the power-law potential 1^(0) = C0^. The 

conditions (jl.9|) and ()1.10|) for the existence of de Sitter fixed points (i^o, 0o) yield 

^°-6(i+e0g)- 12^ ^ ' '^°-(p-4)r ^^-^^ 

These expressions correspond, e.g., to eqs. (3.4) and (3.5) of Ref. [32] • If ^ < it must 
be < p < 4 for the existence of these fixed points while if ^ > it must be p < or 
p > 4. The stability condition ()2.3j) yields, assuming a positive effective gravitational 
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coupling, i.e., 1 — ^ 0q > 



< 7^-^ • (3-2) 



By using the expression of 0o ()3.1|) and remembering the condition for the existence 
of de Sitter space, one concludes that: 



• If ^ > and p < the fixed point {Hq, 0o) is stable. 

• If ^ > and p > 4 this fixed point is stable. 

• If ^ > and < p < 4 this fixed point is unstable. 

• If < and < p < A this fixed point is unstable. 



Therefore, de Sitter attractors exist in phase space for > and for p < and p > 4. 
The requirement that the de Sitter point be an attractor does not constrain ^ apart from 
its sign, but the spectrum of density perturbations generated during inflation does set 
constraints on ^ for various choices of the power p . 



4 Induced gravity inflation 

We now turn our attention to induced gravity which has been used as the context 
for both dark energy and inflationary models with power law potential [321 EZI- This 
theory is described by 

/(0,i?)=e02i?, ^ = 1, V{(P) = Cr- (4.1) 

Conformal coupling corresponds to e = —1/6. Here we restrict ourselves to considering 
positive gravitational coupling e > 0.^ By using eqs. ()1.9p and p.lO|) one finds that 

. (P^^ (4.2) 

and it would appear that any power p in the potential gives an admissible de Sitter 
space, but further comparison with eq. (jl.9|) . which assumes the form —SeH^ +C0q = 
yields the result that de Sitter equilibrium points {Hq, 0o) can only exist if p = 4. In 
this case, a mass term m^cjp' /2 can not be added to the the potential if one wants to 

^For negative gravitational coupling e < de Sitter solutions with the quartic potential V ~ Xcf)^ /A 
exist only if A < 0. 
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find de Sitter solutions and use them as infiationary attractors. For V{(j)) = A0'^/4 the 
de Sitter fixed points are given by^ 




A 02 



12e 



,00 



(4.3) 



including the Minkowski space (i^o^^o) = (0,0). Discarding the contracting de Sitter 
spaces with Hq < that are always unstable [22] , the stability condition p. Ill) reduces 
to 12eHQ — A0Q < 0. By using eq. (j4.2|) it is seen that this stability condition is always 



respect to inhomogeneous perturbations. To the best of our knowledge, only stability 
against homogeneous perturbations was considered in the literature. 

5 Conclusions 

The criterion p.llj) for the stability of de Sitter space in scalar-tensor gravity allows one 
to quickly assess stability of de Sitter fixed points given the scalar-tensor Lagrangian of 
the theory, once conditions ()1.9|) and p.lOj) establish their existence. The advantage of 
this stability condition is that it establishes stability with respect to the more general 
inhomogeneous perturbations and already incorporates a covariant and gauge-invariant 
treatment. By contrast, most studies of stability available in the literature are restricted 
to the less general homogeneous perturbations. The main limitation of the stability 
criterion p.llj) is that its validity is restricted to linear perturbations. While this suffices 
for most situations, second or higher order corrections would be needed, e.g., to analyze 
the backreaction of long wavelength perturbations on the de Sitter metric. 

The simplification introduced by (jl.llj) in the stability analysis of de Sitter space 
is made clear by the scenarios discussed in this paper in which we can quickly recover 
the presence of de Sitter attractors in one scenario of extended quintessence and in two 
scalar-tensor infiationary scenarios, and extend the results of previous studies. There 
are other scenarios in the literature for which one would like to be reassured in a gauge- 
independent way about the stability of a crucial de Sitter attractor with respect to the 
more general inhomogeneous perturbations. These scenarios will be analyzed in future 
pubhcations. 

^This expression matches, e.g., eq. (7d) of Ref. . 



satisfied for e > and A > and the fixed point ( y 0o 




is stable, to first order, with 
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Appendix: Derivation of the stability criterion ( 11.111) 

Here the stability criterion (jl.llj) is derived using the covariant and gauge-invariant 
formalism of Bardeen in the form developed by Ellis, Bruni Hwang, and Vishniac |23j . 
A version of this formalism directly applicable to scalar-tensor gravity was given in 
Refs. [21]. The metric perturbations are defined by the equations 



9oo 



-a^ (1 + 2AY) 



(A.l) 



9oi 



-a'BYi 



(A.2) 



9ij = a' [in (1 + 2ifL) + ^Ht Yij] , 
where the scalar harmonics satisfy the eigenvalue problem 

V,- V* Y = -k^Y , 



(A.3) 



(A.4) 



7ij are the components of the three-dimensional FLRW background metric, and Vj is the 
covariant derivative of 7jj. The defining relations for the vector and tensor harmonics 



Yi and Yij are 



Y 



1 ^ 



1 ^ ^ 



(A.5) 



The variables used are the Bardeen potentials and $yi and the Ellis-Bruni variable 
A$ given by [22] 



k 



a • / a • 



(A.6) 

(A.7) 
(A.8) 
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The gauge- independent variables AF, A/, and AR are introduced by equations similar 
to eq. (rO|l . 

By specifying that the background spacetime is a de Sitter one, the perturbations 
satisfy, to first order, the equations 

A(j) = 6(j), AR = SR, AF = SF, Af = 5f. (A.9) 

The first order equations ruling the evolution of these perturbations reduce to 



A0 + 3HoA(p + 



a? 2ujq 



(/o - 2^,") 



f4>R 
2uJn 



AR 



AF + 3HoAF+ { — -AHn AF + ^AR = 



Hx + ^Hq H — ^7 Hx = 



I ( AF Af' 
+ Ho^A = - \ - Ho — 

^ \ -To -To 



1 AF 



H 



2 Fo 
AF 



"0 



(A.IO) 

(A.ll) 
(A.12) 

(A. 13) 
(A.14) 
(A.15) 



where Ai? is now given by 



AR = 6 



2k^ 



lAF AF 3H^ AF 



3 a 



and with the notations 



d(t)dR 



V3a2 



{4>o,Ro) 



(A.16) 



(A. 17) 



(A.18) 
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By comparing eqs. ()A.14|) and ()A.15|) one obtains 

AF 



<I>H = $A 

Further substitution in eq. ()A.17|1 yields 



2Fn 



AR = 6 



(A. 19) 



(A.20) 



The stabihty with respect to tensor perturbations was studied in Ref. j22j, with the 
result that de Sitter space is always stable, to first order, with respect to these modes. 
The vorticity modes cannot be generated in scalar-tensor (or generalized) gravity theories 
when matter contributions are absent. There remain the scalar modes. We now restrict 
our attention to scalar-tensor gravity, for which / (0, R) is linear in the Ricci curvature 
R. By using ^ = ^ A0, eq. (|Xl9|l gives (see Ref. |22] for further details) 



H 



2 Fo 



A(j) 



while 



AR 



•^f<f>R 



A<p + 3iJo + 



Acf) 



(A.21) 
(A.22) 



from eqs. ()A.20|) and ()A.21|) . This, and eq. ()A.10|) . in turn yield 
A0 + A0 



7^ 



2 "^O + Fo 



6/|« ^2 



^0 1 + 



3/, 



<t>R. 



A0 = 



(A.23) 



when 1 + 3f^J{2iUoFo) 0. When instead 1 + 3/|r/(2cuoFo) = 0, eq. (OrTO|l has either 
the trivial solution A0 = or - 2Vq = Suq 
The choice 



A0 = ee^*, (A.24) 

provides solutions of the asymptotic form of eq. ()A.23|) at late times. Here e and 7 are 
constants and 7 obeys the algebraic equation 



72 + 3i/o7 + c = 



with 



f" 

2 ''^0 



(A.25) 



(A.26) 
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The roots of of eq. ()A.26|) are 



It is straightforward to see that Re (7„) < 0, while the sign of Re (7+) depends on the 
sign of c. If c > then Re (7+) < and we have stability. If instead c < then 
Re (7+) > and de Sitter space is unstable. The condition for the stability of de Sitter 
space in the scalar-tensor theory p.5|) with R) linear in R is therefore 



Fo ^0 



< . (A.28) 



2(^0 Fo J 
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